Quantum tetrahedral mean field theory of the magnetic susceptibility for the pyrochlore lattice 
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Abstract 

A quantum mean field theory of the pyrochlore lattice is presented. The starting point is not the individual magnetic ions, as in the usual 
Curie-Weiss mean field theory, but a set of interacting corner sharing tetrahedra. We check the consistency of the model against magnetic 
susceptibility data, and find a good agreement between the theoretical predictions and the experimental data. Implications of the model and 
future extensions are also discussed. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.40.Cx 



Introduction.- Geometrically frustrated antiferromagnets 
with a pyrochlore lattice exhibit a rich phenomenology which 
has received a vast amount of attention during the last decade 
0]. In the pyrochlore lattice, the magnetic ions occupy 
the corners of a 3D arrangement of corner sharing tetrahedra. 
Materials that crystallize in this structure exhibit anomalous 
magnetic properties ^ ||]: The magnetic susceptibility fol- 
lows the Curie- Weiss law down to temperatures well below 
the Curie temperature. Some of them exhibit long range or- 
der at very low temperatures, whereas others behave as spin 
glasses, even though the lattice is almost perfect. There are 
other compounds of this class that show short range order, and 
are regarded as spin liquids. 

In spite of the intensive research in these systems during 
the last decade, there is still a lack of critical comparison be- 
tween theory and experiment. Even though there are a number 
of models and classical Monte Carlo simulations [|| which 
qualitatively describe some of the experimental results, some 
of the features found in the experimental data cannot be ex- 
plained by means of a classical theory, as are, for example, the 
maxima appearing in the magnetic susceptibility at very low 
temperatures, which a classical model cannot explain. How- 
ever, there have been few attempts to investigate these systems 
quantum mechanically. Harris and co-workers have studied 
the quantum s — ^ Heisenberg antiferromagnet [Q]. Canals 
and Lacroix [|§]], have applied a perturbative approach to the 
density operator of a small cluster, and found that the ground 
state is a quantum spin liquid. 

In this work, we undergo the task of making such a quantum 
theory of the pyrochlore lattice in the framework of the mean 
field theory. The goal of this work is twofold: In one hand, we 
introduce a fully quantum mechanical mean field theory of the 
Heisenberg antiferromagnet in the pyrochlore lattice for arbi- 
trary spin s. In the other hand, we try to make this model as 
simple as possible, in order to make it easy to compare it with 
experimental data and extract information about the various 
interactions that play a role in these systems. Only the main 
results of the model will be presented here, as a more detailed 
presentation of the model will be published elsewhere. The 
starting point to build this MF is not a set of interacting spins, 
but a set of coupled tetrahedra, which goes back to the con- 
stant coupling approximation of Kastelejein and van Kranen- 
donk [0]. The magnetic susceptibility of this system is cal- 
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TABLE I: Values of g(S) for the individual spins, s, considered in 
this work. 



culated for both the interacting and non-interacting tetrahedra 
cases. We find that the independent tetrahedra model fails to 
explain the behavior of experimental data, and we justify this 
fact, whereas the MF theory describes quasi-quantitatively the 
experimental data for a variety of systems. 

The model.- The Hamiltonian of the quantum Heisenberg 
model with nearest neighbors in the presence of an applied 
magnetic field Hq in the pyrochlore lattice can be put as JTo| ] 



H 
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(id) 



(1) 



where J\ is the negative exchange coupling, and Si is the spin 
operator of the i-th magnetic ion. The Zeeman term has been 
quoted in units of the Bohr magneton times the gyromagnetic 
ratio, so Hq has dimensions of energy. 

We start by considering the magnetic susceptibility of one 
tetrahedron. In this simple case, the Hamiltonian can be easily 
diagonalized in terms of the total spin representation of the 
tetrahedron, and the magnetization can be expressed as 



M 



J2 s g(S)Se^ s +VQs(x) 
Ecg(g)e^ s ( s + 1 ) Sinh( 7f^ ' 



(2) 



where S represents the modulus of the total spin operator of 
the tetrahedron; g(S) is the degeneracy associated to the total 
spin value S, which can be calculated by using Van Vleck's 
formula [fill], and are listed in table | for the values of the 
individual spins s considered in this work; x = ^Sf-, and ji = 
ik (the energies are quoted in units of the Boltzmann constant, 



1 



so they have dimension of absolute temperature); and 



Qs(x) 
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In the limit x <C 1, we can define the magnetic susceptibil- 
ity of the tetrahedron as 



M _ 1 J2s g(S)S(S + 1)(25 + l)e^ s ( s+1 ) 



Es9(S)(2S+l)e^(s+i) 



(4) 



The susceptibility per ion in the tetrahedron will be simply 



given by x tet = 2L ^-. 

It is interesting to note that for 
ity can be identified with the high 
a Curie-Weiss type law (per ion 
are considering in this work), x tet 
9 = 2s(s + 1)J\. Therefore, this 
havior predicted by the Curie law 
However, the value of the Curie- 



ji <C 1, the susceptibil- 
temperature expansion of 
and in the same units we 
= £l£+Ii + where 
model reproduces the be- 
at very high temperatures. 
Weiss temperature in this 



type of lattice is given by [ 10 1 



U%s(s + 1) = 29. 



(5) 



The reason for this deviation is, obviously, the fact that we 
have not considered the interaction with the neighboring tetra- 
hedra, which is of the same order of magnitude that the inter- 
actions inside the tetrahedron. 

To remedy this situation, we can introduce a tetrahedral 
mean field (TMF) theory that takes into account, at least in 
an approximate way, the interaction with the neighboring ions 
outside the tetrahedron. Each ion in the tetrahedron interacts 
with z = 3 external ions, in the nearest neighbor interaction 
approximation. 

The effect of this interaction with the nearest neighbors out- 
side the tetrahedron, in the spirit of the MF theory, can be ac- 
counted for by introducing a molecular field proportional to z 
and the magnetization per ion m = The constant of pro- 
portionality can be put without loss of generality as 2J e ff. In 
the high temperature limit, the magnetization per ion can be 
put asm = x* et (Ho + 2zJ e f fin), from which we obtain the 
expression of the susceptibility in this tetrahedral mean field 
(TMF) model 



tmf 



X 



1 - 2zJ eff x 



CAet ' 



(6) 



The value of the J e / / parameter can be estimated as fol- 
lows: in the high temperature limit, x* can be put again 
as the high temperature expansion of a Curie- Weiss type law. 
By equating this expansion with the real Curie- Weiss law up 
to 7^2 terms, we reach at the condition 



9 + 2zJ eff 



8(8 + 1) 



(7) 



Were the interaction with further neighbors negligible with 
respect to the first neighbor interactions, relation (^) would 
be exact and, therefore, by using the relation between 9 and 
J\ obtained in the non-interacting tetrahedra case, we would 
have J e ff = J\, which is completely similar to the value ob- 
tained in the standard mean field theory. 

Of course, second and further nearest neighbor interactions 
are always present in real systems. Thus, the value of 9 
obtained from experimental data can contain additional con- 
tributions coming from next nearest neighbors and so on. For 
this reason, relation (^|) is only an approximate one, and so 
is ([7]). In order to compare the predictions of the model with 
experimental measurements, it is preferable to consider the 
interaction inside the tetrahedron, J\, and the interaction with 
nearest neighbors outside the tetrahedron, J e ff, as adjustable 
parameters. The difference between J\ and J e // provides 
a way of estimating the value of additional interactions not 
explicitly accounted for in this model. In fact, interactions 
with further neighbors can be explicitly accounted for in this 
model in the following way: let us consider the interaction 
with z' next nearest neighbors (9 for the pyrochlore lattice) in 
the mean field approximation. In this case, the magnetization 
per ion will be given by m = x tet (Ho + 2zJ\m + 2z' </ 2 ra), 
where J 2 represents the coupling with next nearest neighbors, 
and the magnetic susceptibility is given by 



X 



tmf _ 



X tet (Ji,T) 



l-6(J 1 + 3J 2 )x tet (JuTy 
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where we have made use of the fact that for the pyrochlore 
lattice z = 3 and z' = 9 and that the difference between 
Ji and J e f f comes from these additional interactions, that is, 

Jeff = J\+ 3.-72- 

Comparison with experimental data.- In order to test the 
validity of the above model, it is important to compare its 
predictions with magnetic susceptibility data. We have se- 
lected the following systems whose magnetic ions form a py- 
rochlore lattice: ZnCr 2 4 [[§, Gd 2 Ti 2 7 [[§, Y 2 Mo 2 7 [jig], 
and CdFe 2 4 @], where the value of the magnetic ion spin 
is s = |, s = |, s = 1, and s = |, respectively. 

The temperature dependence of the magnetic susceptibility 
in these materials can be seen in fig. []]. All of them follow 
a Curie- Weiss law from high temperatures to well below the 
Curie- Weiss temperature. The predictions of our TMF theory 
are also plotted in fig. [j] The values of J\ and J 2 , which we 
obtain by fitting the data points above T c , have been quoted in 
table ||. Here T c marks the transition to the low temperature 
phase, which is either long range order (LRO) or spin-glass 
(SG). The data for Gd 2 Ti 2 C<7 are from AC susceptibility mea- 
surements. Also, the predictions of the non-interacting tetra- 
hedra model have been represented for the case of ZnCr 2 C»4, 
in order to show the deviations of the predictions of this model 
from the experimental data. The case of CdFe 2 04 is some- 
what special and will be analyzed below. 

As we can see from observation of the figure and the values 
collected in table the agreement between theory and exper- 
iment is quite good. In the case of ZnCr 2 04, the agreement 
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FIG. 1: Comparison with experimental susceptibility data of (a) ZnCr204, (b) GcbT^Oy, (c) Y2M02O7, and (d) CdFe204. The open circles 
represent experimental data, whereas the solid lines are fits to expression (^). The dashed lines represent fits to the Curie-Weiss law. The 
dot-dashed line in the case of ZnCr204 represents the prediction of the isolated tetrahedra model with the parameters obtained from the fit to 
the Curie-Weiss law. 



is excellent in the entire temperature range. For the other two 
systems, there is only qualitative agreement. Nevertheless, 
we feel that the agreement between theory and experiment is 
good enough to assess the validity of our model, at least, from 
a qualitative point of view. 

It is important to notice that the values predicted for Ji are 
systematically smaller than the values of J e ft = J\ + 3J2, 
which are close to the ones obtained from fits to a Curie- Weiss 
law, which is consistent with the idea mentioned in the previ- 
ous section about additional interactions from further neigh- 
bors. 

An additional prediction of our model is the existence of 



a maximum appearing in the \ versus T plots. The position 
of this maximum, T m , presented in table [5], can only be cal- 
culated numerically, as its determination involves solving a 
transcendental equation. However, it is found that it follows 
this empirical law 



T„ 



(a + 1) Ji 
1.273 



(9) 



It is important to notice that the position of the maximum is 
determined only by the value of the interaction inside the tetra- 
hedron. We see that the value of T m predicted by the model in 
the case of ZnC^C^ is very close to the experimental value. 
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However, in the case of Gc^l^Oy the predicted maximum 
occurs at ~ 1/2 of the experimental value. In the case of 
Y2M02O7 there is no maximum in the paramagnetic phase as 
there is a transition to a SG at around 18 K. 

As commented above, the case of CdFe2C>4 is especial, be- 
cause in this material, next nearest neighbor interactions are of 
the same order of magnitude than the nearest neighbor ones, 
but ferromagnetic, which leads to an experimental value of the 
Curie- Weiss temperature equal to zero. However, the present 
model is able to describe even this limiting case, and the val- 
ues of the various couplings obtained from the fit are very 
reasonable. Moreover, the position of the maximum predicted 
by the model is in very good agreement with the experimental 
value, and inside the experimental uncertainty. 

Conclusions. - In this work we have analyzed the tempera- 
ture dependence of the magnetic susceptibility data of an ar- 
rangement of magnetic ions of spin s in a pyrochlore lattice. 
To do this, we have developed a tetrahedral mean field the- 
ory taking as a starting point the exact susceptibility of a set 
of four interacting spins in the corners of an isolated tetrahe- 
dron. For the sake of completeness, we have also briefly an- 
alyzed the susceptibility predicted by a non-interacting tetra- 
hedra model. We reach the conclusion that, even though the 
non-interacting model provides an adequate description of the 
high temperature region, it fails to quantitatively describe the 
temperature dependence of the susceptibility at lower temper- 
atures, as the term in this model is | of the actual value. 
However, once we incorporate the interactions with the rest 
of nearest neighbors in terms of a mean field, we find a quite 
good agreement between the theoretical predictions and ex- 
perimental data. Especially important is the prediction of the 
appearance of a maximum in \ versus T. The position of the 
maximum has been calculated, and a reasonable agreement is 
found. At this point, it is very difficult to say if the deviations 
are an intrinsic problem of the model. To clarify this point, 
more good quality experimental data in a variety of systems 
whose susceptibility maxima lie well above T c would be nec- 
essary. 

One feature of this model is that it can be very easily gen- 
eralized to include additional interactions coming from fur- 
ther neighbors and more exotic systems where, for example, 
the interaction with the nearest neighbors is antiferromagnetic 
and is ferromagnetic with next nearest neighbors, as we have 
done for the case of CdFe2C<4, finding a very good agreement 
between theory and experiment. 

Of course, this work is not conclusive, in the sense that it 
does not solve the problem of the anomalous behaviors found 
in these types of systems. It does not provide an explanation of 
why some of these materials exhibit a long range order state at 
very low temperatures, whereas other are in a spin liquid state. 
Moreover, nothing has been said about the disorder always 
present in any material, which has been suggested to be related 
to the appearance of behaviors characteristic of spin glasses. 



These subjects are out of the scope of this paper, though they 
should constitute the direction of future work. 

In any case, we think that the present model provides a good 
starting point for such investigations. 
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Material 


qCW 


jCW a 


rpexp b 
-*- m 


r c 


Ji 


J 2 


rpfit C 


ZnCr 2 4 


-388 


-25.9 


37.4 


12 (LRO) 


-19.2 


-2.3 


38.0 


Gd 2 Ti 2 7 


-21 


-0.33 


2.0 


0.97 (LRO) 


-0.20 


-0.05 


0.7 


Y 2 Mo 2 7 


-61 


-7.625 




18 (SG) 


-7.13 


-0.02 


11.2 


CdFe 2 4 








17.0 


10 (LRO) 


-5.23 


4.1 


14.4 



"Value of the antiferromagnetic exchange coupling obtained from the equa- 
tion 6 CW = 4J cw s(s + 1) 
^Experimental value o£.the position of the maximum 
c Calculated from eq. (H) 

TABLE II: Parameters of the TMF theory. All of them are quoted in K. 
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